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I.  INTRODUCTION 


The  projections  of  antenna  pattern  characteristics  into  contours  of 
constant  power  flux  on  the  earth  are  called  footprints.  The  coverage 
characteristics  of  these  footprints  are  used  in  studying  satellite 
communication  links,  to  infer  the  (roll,  pitch  and  yaw)  variations  in  coverage 
as  satellites  depart  from  their  nominal  attitude,  and  in  remote  sensing 
applications  to  determine  the  ability  to  resolve  surface  features.  The 
footprints  are  not  only  a  function  of  the  antenna  pattern  but  also  depend  In  a 
complicated  way  on  the  relationships  or  mapping  between  a  coordinate  system 
fixed  to  the  antenna  and  the  earth's  coordinate  system.  The  coordinate  system 
on  the  antenna  may  also  be  different  from  the  coordinate  system  used  to 
describe  the  location  of  the  satellite  itself. 

Since  the  kinetic  energy  of  any  satellite  by  definition  is  insufficient 
to  escape  the  gravitational  field  of  the  primary  body,  its  unperturbed  orbit 
is  periodic  and  closed.  The  period  of  revolution  determined  by  Kepler's  third 
law  in  the  reference  frame  is 


T 


2 


Ha3 

u 


(l) 


where  a  Is  the  semi-major  axis  of  the  elliptic  orbit  and  p  equals  the  product 
of  the  mass  of  the  primary  body  and  the  universal  gravitation  constant,  G. 

For  the  earth,  a  currently  accepted  value  is 


p  -  3.98600800...  x  101A  m3/sec2.  (2) 

A  synchronous  satellite  has  a  period,  T,  equal  to  the  sidereal  period  of 
rotation  of  its  primary  body;  i.e.,  the  time  during  which  the  earth  makes  a 
complete  rotation  relative  to  the  vernal  equinox.  Note  that  the  period  for 
the  earth  is  not  24  hours,  because  in  one  day,  the  earth  rotates  around  its 
polar  axis  and  also  completes  1/365.24  of  the  annual  earth  orbit  around  the 
sun.  Consequently,  a  geosynchronous  satellite  has  a  period  given  by 
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86164.09054  ephemeris  seconds 


(3) 


If 


^  23h56™  48 


mean  solar  day 


T  -  a  T  (4) 

P 

where  a  <  1  the  satellite  Is  called  subsynchronous. 

In  addition  to  the  shape  of  the  satellite  antenna  beam,  the  1/R  factor 
for  power  density  dependence  can  modify  the  shape  of  a  footprint  If  the  latter 
corresponds  to  a  constant  power  contour.  However,  Figure  1  shows  that  the 
most  this  variation  can  be  Is  about  1.3  dB  for  the  synchronous  altitude.  That 
is, 


20  log ,_(R  /R  .  )  -  1.3  dB. 
s10  max  min 

At  subsynchronous  altitudes,  the  variation  In  range  becomes  more 
pronounced . 

Another  possible  modification  of  the  footprint  Is  caused  by  bending  of 
the  rays  as  they  pass  through  the  atmosphere.  The  effect  would  be  greatest 
for  stations  on  the  horizon;  i.e.,  stations  which  are  required  to  operate  at 
low  elevation  angles.  For  example,  consider  Anchorage,  Alaska,  at  about  60°N 
latitude.  Figure  2  shows  the  geometry  for  estimating  the  amount  of  bending  of 
the  rays.  Thayer  (1961)  gives  the  following  formula  for  the  amount  of  bending 
of  a  radio  ray  in  radians  as  it  passes  through  an  exponential  atmosphere: 

t(rad.)  ■  0.0003  coto  ,  a  >  15®  (5) 

o  o 

where  o  Is  the  take-off  angle  (electron  angle),  which,  from  the  geometry  in 
o 

Figure  2,  is  about  22°.  This  value  for  the  take-off  angle  gives  a  total 
bending  of  about  7.4  x  10”^  radians  or  about  0.043  degrees.  The  take-off 


angle  should  be  greater  than  about  5°  in  a  satellite-communications  link,  to 
avoid  excessive  bending  of  the  rays. 

A  third  modification  of  the  antenna  footprint  is  due  to  the  oblate 
spheroidal  shape  of  the  earth.  Again,  this  effect  is  on  the  order  of  a  few 
railliradians.  A  fourth  modification  is  the  effect  of  atmospheric  absorption 
which  becomes  significant  at  EHF  frequencies.  For  example,  at  10  GHz  the 
atmospheric  attenuation  will  be  about  2  dB  at  the  horizon  (Bean  and  Dutton, 
1966).  Jacobs  and  Stacy  (1971)  develop  expressions  for  computing  the 

footprints  of  satellite  antennas  with  circularly  symmetric  patterns  which 

_o 

takes  into  account  not  only  the  R  factor  but  also  atmospheric  attenuation. 

The  footprint  defines  the  area  on  the  surface  of  the  earth  inside  of 
which  the  system  has  a  specified  or  greater  level  of  sensitivity.  The 
mathematical  problem  of  finding  the  locus  of  the  footprints  can  be  solved  in  a 
straightforward  manner  using  vector  analysis.  This  report  documents  a 
computer  program  written  to  draw  the  locus  of  the  footprints  on  computer 
generated  maps.  A  particular  feature  of  this  analysis  is  the  capability  to 
generate  footprints  for  a  satellite  which  uses  repetitively  scanned  beams  to 
achieve  the  total  required  coverage. 


II.  VECTOR  ANALYSIS 


A.  SUN- SYNCHRONOUS  ORBIT,  ANTENNA  SCANNING  WITH  RESPECT  TO  BODY- FIXED, 

VEHICLE-CENTERED  COORDINATE  SYSTEM 

In  this  section,  the  equations  for  the  locus  of  intersection  points, 
i.e.,  the  footprint  of  the  satellite  antenna  beam  and  the  spherical  earth,  are 
derived.  The  particular  orbit  investigated  is  a  “sun-synchronous"  orbit  for 
which  the  satellite  orbit  will  pass  over  a  given  latitude  at  the  same  local 
time-of-day  for  each  pass.  Because  of  this  similarity  of  lighting  and  the 
near-consistency  of  the  local  time-of-day,  most  meteorological  satellites  like 
the  Defense  Meteorological  Satellite  Program  (Hollinger  and  Lo,  1984)  are  in 
this  type  of  orbit.  The  inclination  of  the  orbit  with  respect  to  the  equator 
in  the  latitude-longitude  coordinate  system  is  shown  in  Figure  3.  In  the 
example  considered  in  this  report,  i  -  98.1°  and  the  circular  altitude  above 
the  equator  is  about  833  km  (Denner,  1982).  The  orbital  period  is  about  101.5 
minutes  and  the  circular  orbital  velocity  is  about  6.58  km/sec,  and  the 
"grazing"  orbit  velocity  is  about  7.406  km/ sec. 

The  geometry  of  the  problem  is  displayed  in  Figure  3.  In  the  derivation, 
the  earth  has  unit  radius  so  all  distances  are  in  earth  radii.  The  position 
of  the  satellite  in  subsynchronous  orbit  from  Figure  3  is  given  by  the 
vector  S, 


S.  =  s  (e^sinScosiJ)  —  e^sin0sini}>  +  e^cos8)  (7) 

where  unit  vectors  are  denoted  by  e  throughout  this  report  and  with  s 
determined  from  the  Vis-viva  ("living  force",  Escobal ,  1965)  equation  for  a 
circular  orbit 


v2  =  u/s  (8) 

with  y  «  1.000952348,  and  v  the  circular  velocity,  and  s  in  earth  radii  units 
(1  circular  satellite  unit  «  7.905293  km/sec.).  From  Figure  3, 
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I, 


e  •  e  ■  cosH  *  sin0cos$ 


(e  x  e  )  •  e  «  sinftsini  ■  cose 
— s  — x  — y 


(e  x  e  )  •  -e  *  sinftcosi  -  sin6sin<|)  (11 

— s  — x  — z 

Basically,  s,  i  and  time,  t,  in  Figure  3  determine  the  orbit  with  the 
angular  variable,  fl,  given  by 

*  <ut  *  vt/s  *  /p  t/s3/2  (12 

From  Figure  3,  a  vector  from  the  satellite  to  an  arbitrary  point  on  the  scan 
track  is 


R  -  e  -  S 

_  _p  _ 


with  the  latitude  (y  -  0q)  and  longitude,  of  a  point  on  the  earth  found 


e  *  e  sin0  cos8  +  e  sin0  sin*  +  e  cos0  (14) 

— p  —x  o  o  — y  o  o  — z  o 

The  vector  R,  in  terras  of  a  coordinate  system  centered  on  the  satellite  with 

the  unit  vector  eT  pointing  in  the  direction  of  motion  and  e  chosen  to  form  a 
1  p 

triad,  is 


R  *  D  f-cosoe  +  sino(sin0e  +  cosf3  e_)  1 
—  1  — s  — p  — T 


cm 


VV.V-. 


(15) 


with  8  -  2it  radians/T,  where  T  is  the  period  to  make  a  scan,  and  o  the  angle 
between  the  boresight  direction  of  the  antenna  and  a  radial  vector  from  the 
center  of  the  earth  and  the  satellite.  From  Equations  (14)  and  (15)  and  the 
fact  that  e  •  e  *  1  we  find 

— p  — p 

no  2  2  2  2 

1  *  [cos  o  +  sin  o(sin  8  +  cos  8) 1  +  2s  D  R  »  e^  +  s 


2  2 
D  -  2sD  coso  +  s 


-  1 


(16) 


and  solving  for  D 


2  2  1/2 

D  ■  s  coso  -  [l  -  s  sin  a]  (17) 

the  negative  root  selected  corresponding  to  the  point  on  the  earth  closest  to 
the  satellite  or  on  the  "lit"  side.  For  example,  if  o  *  ir/4, 


and  if  s  -  1.131463449  e.r. 


D 


s  -  /2-s2 

/2 


and  from  Figure  4 


D  -  0.200 


2  2 

D  ■  1  +  s  -  2s  cos0 


(18) 


(19) 


or 


cos0  „  8  ±  1  -  -P-  -  0.989957 

2s 


6  -  8.127° 


(20) 
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To  find  the  unit  vector,  ^T,  in  the  direction  of  motion, 
Oxyz  about  the  x-axis  by  the  angle  as 


x' 

y' 

z' 


o 

o 

1 _ 

x 

0  sini  cosi 

y 

0  -cosi  sini 

z 

and  the  unit  tangent  vector,  _ej,  in  the  xyz  frame  is 


e_  -  -  e  sinQ  +  e  cosft 
— T  —x  — z 


Therefore  ej  in  frame  x'y'z'  is  from  Equations  (21)  and  (22) 


1 

0 

0 

■  — 

-e  sinfl 

—X 

-T  " 

0 

sini 

cosi 

0 

0 

-cosi 

sini 

e  cosft 
—z 

e_  *  -  e  sinft  +  e  cosfl  cosi  +  e  cosjlsini 

— T  —X  — y  — z 


Then 


,  the  third  vector  in  the  triad,  e^,  e^  is  defined  by 


e„  *  e_  x  e 
— 6  — T  — s 


and  substituting  Equations  (7)  and  (23)  into  (24)  gives 


e„  -  e  sini  -  e  cosi 

-6  -y  -z 


A  check  on  the  unit  vectors  _ej  and  e^  is  provided  by  Equations 
which  implies 


rotate  frame 


(21) 

(22) 


(23) 

(24) 

(25) 

(13)  and  (16) 
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|R  +  s| 


D  +  s  -  2Ds  cosa 


( 


Substituting  for  _R  from  Equation  (15)  and  from  (23)  and  (25),  in  terras  of  ey 
e,  and  e_  shows  indeed  Equation  (26)  is  satisfied. 

Therefore,  we  can  find  typical  colatitudes,  9q,  and  longitudes,  alo 
the  scan  track  from 


and 


S  +  R 


e 

“P 


(2 


2  2 

/D  +  s  2Dscoso 


cos  9  m  e  •  e 
o  — p  — z 


sin  8  sin*  -  e  •  e 
o  o  — p  — y 


sin  0  cosA  *  e  •  e 
o  o  — p  —x 


(2 


Knowing  0q 
footprint ; 
defined  as 


and  $  we  can  now  construct  a 
o 

i.e.,  e,  e  and  e  in  Figure 


rectangular  frame  for  defining  the 
3.  The  first  unit  vector  e  is 


e^  -  R/|R|  (2 

and  from  Equations  (7),  (15),  (23)  and  (25)  and  a  great  deal  of  algebra  one 
finds 


e,  *  -e  (coso  cosft  +  sino  sinfl  cosg)  +  e  (-coso  sinft  cosi  +  sino 
-Z  -x  y 


sing  sini  +  sino  cosg  cosJl  cosi) 


+  e  (-coso  sinfl  sini  -  sino  sing  cosi  +  sino  cosg  cosil  sini) 


From  Figure  3 


e  *  e  coso  +  e  sino 
-n  -5  -s 


(31) 


— 1 —  {-e  f coso  cosfi  +  sino  sinn  cosB  -  — +  e  f-coso  sinn  cosi 
tano  1  — xv  coso '  -vv 


+  sino  sing  sini  +  sino  cosg  cosn  cosi  + 


-y 

sinn  cosi 


coso 


■) 


+  e  (-coso  sinn  sini  -  sino  sing  cosi  +  sino  cosg  cosn  sini  + 


slnn  sini 


cos  o 


)} 


and  the  third  vector  in  the  triad  e^ ,  e  ,  e  is 

"M  — 


e 

-c 


e,  *  e 

-€  -n 


(32) 


e  sing  sinn  +  e  (cosg  sini  -  sing  cosi  cosn) 

-*  -y 


-  e  (cosg  cosi  +  sing  sini  cosn) 
— z 


A  typical  generator  of  the  antenna  footprint  is 


T  ■  t(w)  [e.  +  afe  coso)  +  e  sinu>)] 

-  L-S  “C 


(33) 


where  a  is  the  half-power  antenna  beam  width.  At  a  typical  intersection  point 


which  becomes  after  a  great  deal  of  algebra 


2  .  2st(cosg  +  a  cosusing)  (s  -1) 

c  2  +  2 
1  +  *  (1  + 


(35) 


The  solution  for  t  closest  to  the  satellite  is 


t(u>) 


~-r  -  (-s  (-cosa  +  acosui  sino) 


(1  +  a‘) 


-[s^(-cosa  +  acosa)  sino)^  -  (s^  -  1)(1  +  a^)]^^} 


(36) 


The  colatitude,  0,  and  longitude,  $,  of  a  typical  intersection  point  are 


COS0 


sin0  cosa  -  e  •  e 
T  — p  — x 


sin0  sin a  ■  e  •  e 
-p  -y 


(37) 


and  this  completes  the  derivation  for  the  subsynchronous  case. 

B.  GEOSTATIONARY  ORBIT 

In  this  section  the  equations  for  the  locus  of  intersection  points  of  the 

satellite  in  geostationary  orbit  and  a  spherical  earth  are  derived.  The 

satellite  is  located  in  the  equatorial  plane  (x-y  plane  in  Figure  5)  at  an 

arbitrary  longitude,  X  ,  and  at  an  arbitrary  distance,  s,  from  the  origin.  In 

o 

vector  notation,  the  satellite  position  is 

S  ■  s(e  cosX  +  e  sinX  ).  (38) 

—  — ^  0  — y  Q  ' 

The  maximum  of  the  antenna  spot  beam  points  to  an  arbitrary  longitude  and 
latitude  (X,  $) ,  designated  the  aim  point.  In  vector  notation,  the  aim  point 
is 


19 


e 

-P 


e^  cosXcos$  +  e^  sinXcos$  +  e^  sin^, 


(39) 


which  is  a  unit  vector  since  the  sphere  has  unit  radius. 


In  order  to  relate  the  shape  of  the  antenna  beam  on  the  satellite  to  the 
locus  of  intersection  points  on  the  sphere,  an  orthogonal  coordinate 
system  (e^t  e  .  _e  )  is  erected  at  the  satellite.  From  Figure  5,  a  unit  vector 
from  the  satellite  to  the  aim  point  is  given  by: 


e  -  S 

-5  *  |L  -  S 


(AO) 


Substituting  Equations  (38)  and  (39)  into  (40)  gives 


e_  ■  fe  (cosX  cos*  -  s  cosX  )  +  e  (slnX  cos*  -  s  sinX  )  +  e  sin*]/d. 
— £  x  Y  o  — y  *  o  — z  T  1 


(41) 


with 

7  1/2 

di  “  [l  +  ®  “  2s  cos$  cos(X  -  Xo)]  .  (42) 

When  the  aim  point  and  subsatellite  point  (N  in  Figure  5)  coalesce  so 

that  *  ■  0  and  X  -  X  , 

T  o 

e _($  •  0,  X  ■  X  )  ■  e  cosX  -  e  sinX  (43) 

— C  o  -x  o  -y  o 


since 


dL  -  s  -  1  (44) 

in  thi6  case.  A  second  unit  vector  in  the  right-handed  system  used  at  the 

satellite  is  e  .  Let  e  lie  in  the  r-z  plane. 

— n  — n 


Figure  5.  Geometry  for  the  derivation  of  footprints.  The  geostationary 
satellite  at  S  is  in  the  x-y  (equatorial)  plane. 


e  ■  -e  sin$(cosX  cos$  -  s  cosX  )/d,d- 

— T1  X  O  i  L 

-e^sin$(sinX  cos$  -  s  sinX^/d^d^ 

+e  (d„/d  ) • 

— z  2  1 

In  the  case  $  -  0 , 

«“£-(♦“  0)  » 

— n  z 

as  of  course  is  expected. 

The  third  member  of  the  triad  is  constructed  to  form  a  right-handed 


system;  i.e.. 


Having  taken  care  in  deriving  the  orthogonal  coordinate  system  at  the 
satellite  will  permit  writing  an  expression  for  a  typical  generator  in  the 
beam.  From  Figures  5  and  6 ,  we  can  write 


T  =  t(io)[e^  +  e  (a  cosoj  cosg  -  b  sinai  sing) 

—  -£  — n 

+  e^(a  cosm  sing  +  b  sinm  cosg)],  (52) 

where  t(w)  is  the  length  of  the  generator  of  the  elliptic  cone,  a  and  b  are 
the  semi -minor  and  serai-major  axes  of  the  cross  section  of  the  elliptic  cone, 
respectively,  g  is  the  tilt  of  the  semi-major  axis  of  the  beam  with  respect  to 
the  n-£  coordinate  system,  and  u>  is  the  parameter  that  generates  the  elliptic 
cone  as  u  varies  from  0  to  2 v  radians.  An  elliptic  beam  can  be  generated  by 
using  a  reflector  antenna  with  different  dimensions  in  the  two  principal 
planes( sometimes  referred  to  as  "cut"  parabolas).  The  tilt  is  produced  by 
rotating  the  reflector  about  its  generator  axis. 

From  Figure  5 ,  at  a  typical  intersection  point,  Q,  we  have 

R  -  S  +  T,  (53) 

and  because  _R  lies  on  a  sphere  of  unit  radius, 

R  *  1  =  |JL|2  ■  1  ”  T2  +  S2  +  2T  •  S.  (54) 

Substituting  Equations  (38)  and  (53)  into  Equation  (34)  gives  the  following 
quadratic  equation  in  ”t”  for  the  intersection  points  of  the  antenna  beam  and 
the  earth;  i.e., 
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0 


(55) 


a  -  2s{cosXot(F1/d1)  -  (^V^lV  +  (F3F6/d2)] 
+  sinXQ[I3/d1)  -  (F^/d^)  -  (F^/dj)]} 

2  i 

a2  -  s  -  1 

F^  »  cosX  cos<$>  -  s  cosXq 
F2  -  sin<f>  F ^ 

F^  *  sinX  cos<j>  -  s  sinAQ 


(56b) 

(56c) 

(56d) 

(56e) 

(56f  ) 
(56g) 
(  56h) 

(561) 
(56j  ) 


In  the  computer  program  in  Appendix  A  (F00T2),  the  parameter  Fy  is  tested 
as  w  varies  from  0  to  2*  to  see  if  F?  is  less  than  zero,  since  these  values 
of  u)  correspond  to  antenna  rays  which  do  not  intersect  the  earth.  When  F^ 
equals  zero,  the  antenna  rays  graze  the  earth. 

The  solution  of  the  quadratic  equation  (55)  has  two  roots  corresponding 
to  an  intersection  point  on  the  near  or  lit  side  of  the  sphere  and  an 
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intersection  point  on  the  back  side  or  shaded  side  of  the  sphere.  Since  a^<0, 
the  smaller  of  the  two  roots  or  the  intersection  point  on  the  lit  side  is 


t 


(57) 


and  when  Fy  =  0  or 


a 


1 


- 


(58) 


the  generators  of  the  antenna  beam  are  parallel  to  the  sphere.  Equation  (58) 
defines  the  optical  horizon  or  "limb  line"  to  which  we  return  shortly.  As  an 
example  of  the  solution  of  Equation  (58),  suppose  $  =  0,  \  -  Xq  =  0;  i.e.,  the 
aim  point  is  at  the  intersection  of  the  Greenwich  meridian  and  the  equator. 

The  satellite  is  also  at  the  Greenwich  meridian.  Also,  assume  the  antenna 
beam  is  circular  (N.B.,  in  this  case  the  tilt  angle,  3,  is  superfluous).  Then 


a0  =  1  +  a^, 

*\  *  -2s, 

~  s  ~  l  * 


Then  Equation  (58)  gives 


(60) 


and  for  the  antenna  beam  to  just  graze  the  sphere, 

1 

a  =  /s  -  1  . 

Substituting  (59)  and  (60)  into  (52)  gives 

2  , 

T  =  — - -  [e^  +  -  (e^cosai  +  e^sinw)]  ,  (61) 

/s2  -  1 

the  generator  of  the  limb  line  or  optical  horizon.  Figure  7  shows  the 
distances  comprising  the  grazing  ray  in  Equation  (62).  Note  that  Equation 
(61)  gives 

III  “  /I  •  i  *  i's2  -  1,  (62) 


which  agrees  with  Figure  7. 

The  latitude  and  longitude  for  an  arbitrary  point  on  the  limb  line  for  an 
arbitrary  satellite  longitude,  Xq,  are 

♦  ■  sin  Vs2  “  1  cosm/s)  ,  (63) 


and 


X  “  tan 


-1 


/  2  ,  - 

s  sinX  + -  (-Vs  -  1  sinX  +  slnucosX  ) 

os'  o  oJ 


(64) 


/s2  —  i  2 

s  cosX - f/s  -  1  cosX  +  sinusinX  ) 

os'  o  oJ 


respectively. 


Equation  (52)  represents  the  generator  of  the  footprint,  and  one  can  now 
calculate  the  latitude  and  longitude  of  a  typical  intersection  point  on  the 
footpoint.  From  Figure  5,  we  see  that  the  latitude  of  a  typical  intersection 
point,  is  determined  by 


(N.B.,  _r  is  a  unit  vector).  Substituting  Equations  (38),  (52),  and  (53)  into 
(65)  gives 

4 

-  sin'1[(t(u))/d1)(sinb  +  F5d2)].  (66) 

In  order  to  determine  the  proper  quadrant  for  the  longitude  of  a  typical 
intersection  point,  X^,  one  must  compute 

Fq  =  e  •  R  =*  cos<j>  sinX  (67) 

o  — y  —  i 

and 

Fg  -  •  R  «  cos^  cosXj^  (68) 

Substituting  Equations  (37),  (52),  and  (53)  into  (68)  and  (69)  gives 

F8  -  s  sinXQ  +  t(oi)[(F3/d1)  -  (F4F5/d1d2)  -  (F^/d^l  (69) 

and 

Fg  *  s  cosXQ  +  t(u)[(F1/d1)  -  (^2F5/dld2)  +  (F3F6/d2)]’  (70) 

so  that  the  longitude  is 

=  tan"l(Fg/Fg).  (71) 

In  Appendix  A  the  computer  program  is  given  for  computing  the  latitude 
and  longitudes  (4^,  X^)  of  the  locus  of  intersection  points  at  half-degree 
increments  in  the  parameter  w. 

C.  METHOD  FOR  DETERMINING  WHETHER  ARBITRARY  POINT  ON  SURFACE  OF  A  SPHERE  IS 
INSIDE  OR  OUTSIDE  A  GIVEN  FOOTPRINT 

The  arbitrary  point  may  represent  the  location  of  a  receiving  antenna  on 
the  surface  of  the  earth  and  be  specified  in  terms  of  its  latitude  and 
longitude.  Using  the  "winding  number"  concept  (Ahlfors,  1966),  we  have 


i2ir,  fj  inside 

(72) 

0,  fj  outside 

where  Au  =  ir/360  and  arg  is  the  argument  of  the  difference  of  the  complex 
functions  inside  the  brackets  and 

f(<d)  =  X(tu)  +  i$(uO  (73) 

where  X  and  <j)  are  the  longitude  and  latitude  computed  from  Equations  (71)  and 

(66). 

D.  PROJECTIONS 

The  continental  borders  and  antenna  footprints  are  plotted  according  to 
the  type  of  projection.  That  is,  points  on  the  earth  defined  by  latitude  and 
longitude  are  transformed  to  points  in  the  plane  of  projection  (the  U,  V 
plane)  In  which  the  map  and  footprints  are  plotted.  The  particular  algorithm 
used  to  plot  the  world  map  backgrounds  and  overlay ' data,  such  as  a  satellite 
ground  trace  or  antenna  footprint,  on  the  requested  map  projection  is  WORLD. 

1.  AZIMUTHAL  PROJECTIONS 

The  U,V  plane  Is  tangential  to  the  earth  at  the  point  $  (POLAT,  POLONG), 
which  transforms  to  the  origin,  $'(0,0),  in  the  U,V  plane  In  Figure  8.  Let  P 
be  a  point  on  the  earth  at  an  angular  distance,  A,  from  the  point  <J>  (POLAT, 
POLONG).  Let  B  be  the  angle  between  the  great  circle  $P  and  the  meridian 
at  9.  Then  from  Figure  8 

e_  =  e  cosX  +  e  sinX  (74) 

— R  ^  —x  o  —  y  o 

and 

en  3  e  cosdcosX  +  e  cosSsinX  +  e  sirnS  (75) 

— R^  —x  — y  — z 
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so  that 


with 


T  - 


1^2  -*ll 


COS^ 


\\  -  %  I  -  2  sin  | 


Substituting  (77)  into  (76)  gives 

T 


sinA 


(76) 


Then 


cosB  *  e„  e  ■  sin6/sinA 
— T  — z 


and 


cosA  -  cos5cos(A-X  ) 
o 


with 


R 


.  A 

r  sin  y 


Stereographic  pro  lection 

R  -  TAN( A/2)  -  ( 1  .-COSA) / SINA  As  A  +  180°,  R  -*•  •.  Thus,  the  entire 
surface  of  the  globe  transforms  to  the  entire  U,V  plane.  In  practice, 
distortion  becomes  great  beyond  R*2,  or  A  -  130°. 

Orthographic  projection 

R  -  SINA.  This  projection  plots  a  hemisphere  within  radius  R-l.  I 
maximum  possible  value  of  A»90°. 


Lambert  equal  area  projection 

R  is  calculated  by  the  two  Fortran  statements 

R  -  (l.+COSA)/SINA 

R  =■  2./SQRT(l.  +  R*R) . 

As  A  +  180° ,  R  >  2,  and  the  entire  surface  of  the  globe  is  plotted  within 
a  radius  R-2.  The  maximum  value  of  A- 180°. 

Gnomonlc  projection 

R  *  SINA/COSA  As  A  +  90°,  R  *  “. .  A  hemisphere  is  plotted  over  the 

entire  U,V.  In  practice,  distortion  becomes  great  beyond  R=*2,  or  A  ~  65°. 

Azimuthal  equidistant  projections 

R  *  A  (in  radians)  *  ACOS(COSA)  As  A  180°,  R  ♦  v..  The  entire  globe 
surface  is  plotted  within  a  radius  R  *  it. 

Cylindrical  projections 

The  U,V  plane  must  be  Imagined  to  be  wrapped  around  the  globe  to  form  a 
cylinder,  the  U-axis  touching  the  globe  on  some  great  circle  (see  Figure  9). 
The  axis  of  the  projection  Is  perpendicular  to  this  great  circle  and  parallel 
to  the  V-axis.  The  point  $(P0LAT,P0L0NG)  transforms  to  the  origin,  $'(0,0), 
of  the  projection  lines  on  the  great  circle.  The  limits  of  the  U-axis  are 
defined  by  a  cut  in  the  cylinder  parallel  to  its  axis  and  diametrically 
opposite  to  $.  The  pole  of  the  projection,  Q,  is  the  point  90°  from  the  great 
circle  in  the  direction  of  +V.  ROT  is  the  angle  between  the  V-axis  and  north 
at  These  points  and  N,  the  north  pole,  are  shown  in  Figure  9.  Points  on 

the  surface  of  the  globe  are  transformed  to  points  in  the  U,V  cylinder  by  the 
rule  appropriate  to  the  projection. 

The  latitude  and  longitude  of  Q  are  calculated  in  terras  of  the  latitude 
and  longitude  of  $  and  the  angle  R$T.  The  angle  R$T1  (see  Figure  9)  is  also 
computed . 

In  Figure  9,  P  is  some  general  point  on  the  surface  of  the  globe.  A  is 
the  angular  distance  of  P  from  Q.  B  is  the  angle  between  the  great  circles  QN 
and  QP.  The  quantities 
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SINA  - 

SIN(A), 

C0SA  =  C0S(A) 

SINB  =* 

SINB(B), 

COSB  =  C0S(B) 

SINR  = 

SIN(ROTl) , 

C0SR  =■  COS(ROTl) 

are  computed. 

Let  P'(U,V)  be  the  point  on  the  U,V  plane  corresponding  to  the  point  P  on 
the  surface  of  the  globe.  Then  U  is  proportional  to  the  angle  a  (see  Figure 
9).  The  value  of  V  depends  upon  the  type  of  projection.  The  coordinates  of 


P'  are  given  below. 

Cylindrical  equidistant  projection  U  -  a  (in  degrees) 

»  ATAN2(SIN(B+R0T1),  -  COS( B+R0T1 ) ) /F 

*  ATAN2(SINB*C0SR+C0SB*SINR,SINB*SINB*SINR-C0SB*C0SR)/F 

Here  division  by  F  converts  radians  to  degrees. 

V  *  90. -A  (in  degrees) 

*  90.-ACOS(COSA)/F 

The  entire  surface  of  the  globe  is  transformed  to  a  rectangle  in  the  U,V  plane 

-180.  <_  180. 

-  90.  <  90. 


Mercator  projection  with  arbitrary 


U  ■  a  (in  radians)  _ 

=  ATAN2(SINB*C0ST+C0SB*SINR,SINB*SINR-C0SB*C0SR) 
V  »  AL0G( C0T(A/2)) 

*  AL0G( ( 1 .+C0SA) /SINA) 


The  entire  surface  of  the  globe  is  transformed  to  an  infinite  rectangle 

in  the  U,V  plane.  As  A  +  0,  V  +  »;  as  A  ♦  180°,  V  +  -  «.  When  a  =  180°,  U 

*  it,  when  a  =  -180°,  U  ■  tt  ,  when  a  =  -180°,  U  =  tt.  Hence 
-*  <  V  <  ® 


— TT  <  U  <  7T 
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In  practice  distortion  becomes  great  for  A  <  5°  or  >175 


Mollweide-type  projection 

The  projection  used  is  not  a  true  Mollweide.  The_  coordinates  of  P'  are 
given  by 

V  -  COSA 

U  is  given  by  the  two  Fortran  statements 
U  -  ATAN2( SINB*COSR+COSB*SINR, SINB*SINR-COSB*COSR) 

U  +  U*.5*SQRT(1-V*V)/ATAM1.0) 


The  entire  surface  of  the  globe  transforms  to  an  ellipse  in  the  U,V 
plane.  The  major  and  minor  axes  of  the  ellipse  are  along  the  U  and  V  axes, 
respectively. 

-2.  <_  U  <_  2. 

-1 .  8  V  <  1 . 


III.  EXAMPLES 


The  following  examples  will  serve  to  demonstrate  the  use  of  the  computer 
programs  given  in  Appendix  A  and  provide  insight  into  the  shape  of  footprints 
produced  by  various  beams  on  satellites  in  sunsynchronous  or  geostationary 
orbit.  The  world  map  was  produced  using  subroutine  WORLDS  (Gurlitz,  1981). 

Consider  the  following  example  for  a  sunsynchronous  satellite  located  a 
distance  s  =  1.131463449  e.r.  and  at  an  inclination  of  98.1°  (this  corresponds 
to  a  DMSP  application).  The  half-power  beamwidth  is  0.005069  rad. 
and  o  »  n/4.  Figure  10  shows  the  footprints  at  five  instants  of  time 
separated  125  millisec  (the  entire  scan  in  this  example  takes  1.9  sec.)  and 
for  eight  passes.  Note  how  the  given  circular  antenna  beam  becomes  an 
elliptical  footprint  with  major  axis  nearly  parallel  to  a  meridian  at  t  *  0 
while  the  footprint  becomes  an  ellipse  with  major  axis  nearly  parallel  to  the 
equator  540  millsec.  later. 

Figure  11  shows  the  received  voltage  versus  time  for  the  first  seven 
passes  of  the  sunsynchronous  satellite.  The  period  for  a  pass  is 
1900  millisec.  The  flattening  at  the  top  of  the  beam  In  Figure  11  results 
from  the  limited  sampling  of  the  antenna  pattern  data  near  the  beam  peak.  The 
location  of  the  transmitter  and  2  and  10  dB  footprints  for  the  first  three 
passes  and  part  of  pass  4  are  shown  in  Figure  12.  The  data  for  the  satellite 
antenna  beamwidth  were  taken  from  the  LFMR  Midterm  Status  Review,  1  February 
1985. 

Next,  consider  a  satellite  in  geostationary  orbit  at  6.619  e.r.  and  \  = 

o 

100°W.  The  3-dB  beamwidths  in  the  two  planes  are  0.017454  rad.  and  0.03491 
rad.  respectively  and  the  tilt  angle  8  =  35°.  Figure  13  shows  the  footprint. 
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Figure  10. 


Footprints  for  satelite  in  sun-synchronous  orbit,  at  an 
inclination  of  98.1°.  The  separation  in  time  of 
footprints  in  125  millisec.  The  second  track  occurs 
1900  millisec  after  track  number  1. 
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.  Time  series  for  sun-synchrono-.^  orbit  at  an  inclination 
of  98.1°.  The  ordinate  represents  the  received  signal 
power  from  a  transmitter  or  the  ground  as  shown  in 
Figure  12. 


IV .  CONCLUDING  REMARKS 


The  vector  analysis  plus  numerical  method  for  generating  earth  coverage 
footprints  for  subsynchronous  and  geostationary  satellites  is  given.  The 
computer  algorithms  should  provide  a  useful  tool  for  determining  the  coverage 
of  regions  on  the  earth  by  various  types  of  antenna  beams  and  configurations 
of  antennas. 

These  earth  coverage  footprints  are  used  in  the  study  of  the  requirements 
of  the  spacecraft  antennas.  For  example,  some  of  the  considerations  are:  1)  a 
narrower  spot  beam  on  the  earth  requires  a  larger  satellite  antenna  reflector; 
2)  the  overlap  of  spot  beams  may  result  in  the  need  for  multiple  frequencies 
or  the  use  of  orthogonal  polarizations;  3)  larger  spacecraft  antennas  mean  an 
increase  in  gain  which  in  turn  may  result  in  a  decrease  in  transmitter  power; 
4)  the  beam  coverage  influences  the  number  of  transponders  used  at  the 
satellite;  and  5)  communication  and  radar  systems  and  better  ground  resolution 
for  radiometric  remote  sensing  applications. 

The  question  of  satellite  coverage  on  the  surface  of  the  earth  is  not  new 
(Siocos,  1973)  and  the  derivations  of  the  geometrical  distances  are  probably 
numerous.  However,  the  derivations  given  in  this  report  based  upon  vector 
analysis  are  simpler  than  those  based  upon  spherical  trignometry  (Jacobs  and 
Stacey,  1971;  Adamy,  1974).  Also,  the  versatility  of  the  shape  and  tilt  of 
the  satellite  antenna  beam  (i.e.,  elliptic  in  cross  section  with  arbitrary 
orientation  of  semi-minor  and  semi-major  axes)  is  an  added  feature  not  found 
elsewhere  as  far  as  is  known.  Jacobs  and  Stacy  (1971)  and  Adamy  (1974) 
develop  mathematical  expressions  for  computing  earth  footprints  assuming 
circular  symmetrical  antenna  beams. 
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APPENDIX  A 


COMPUTER  PROGRAMS 


► 

w 

\ 

N 


The  following  is  a  listing  of  the  control  "cards"  for  compiling,  linking 
and  executing  the  program  for  drawing  footprints  for  a  sun- synchronous  satel¬ 
lite.  The  job  is  submitted  as  a  batch  job  to  the  CDC  176  (Mainframe  X). 


100  F19EZ ,STMFX(mainf rame) ,P3000( priority )  ,T5000( time  limit  in  octal  sec.), 
MS160000( storage) ,NT0(no  magnetic  tapes) 

110  ACCOUNT  (OTT,  R.  17533  F19E2574910  A21259  54115  6530) 

120  XMIT, OUTPUT.  (send  output  file  to  IBM  printer) 

130  FILE, FOOT,  TR-Z,BT-C,FL-80 . 

140  ATTACH, FOOT, ID- 175 33, ST-PF6  (attach  FORTRAN  prog) 


150  COPY, FOOT, IN. 

160  REWIND,  IN. 

170  FILE, DATA, RT-Z,BT-C,FL»80. 

1 80  ATTACH , DATA , DATA , ID- 1 7  5  33 , ST-PF6 . 
190  FTN5 , I- IN , LO-R/ A/M/ S , PL-1 5000 . 

200  ATTACH .WRLLIB.WORLD5LIB , ID-12730 . 
210  ATTACH, PLIB.FTN5PLOTLIB. 

220  ATTACH, TAPE8,1BKGD, ID-1 27 30 

230  LIBRARY, WRLLIB.PLIB. 

240  LOAD(LGO) 

250  EXECUTE( .DATA) 

260  HARDCPY . 


FTN5  (FORTRAN'77)  does  not  interpret 
Record  type  -  Z  format.  Therefore, 
need  the  copy  command  which  converts 
FOOT  with  record  type  Z  to  record  type 
W  (word  control).  The  default  record 
type  on  upper  Cyber  is  W. 


Attach  data  file  giving  latitude  and 
longitude  of  central  point  of  projec¬ 
tion,  projection  type,  flag  for  plott¬ 
ing  data  from  Tape  12,  latitude  and 
longitude  grid  spacings,  etc. 


Attach  map  plotting  routine 

:Attach  system  plotting  library 

Attach  continential  background  data 
file 

Link 

LOAD  sequence,  here  the  load  consists 
of  modules  from  three  different  files 
and  execute  as  a  single  program 


Command  for  making  hardcopy  plots. 
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w  ■r. v 


1  oo 

1 10 
120 
130 
140 
130 
160  C 
170  C 
180  C 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
430 
460 
470 
480 
490 
300 
510 
520 
530 
340 
530 
560 
570 
380 
390 
600 
610 
620 
630 
640 
630 
660 
670 
680 
690 


PROGRAM  F00T4 ( INPUT , TAPE5= I NPUT , OUTPUT , TAFE6-QUT PUT ) 
DIMENSION  UFOOT( 181.84) , VFOO T < 1 81 . 84 ) 

DIMENSION  A (30) , ID (7) 

COMPLEX  W  < 1 81 ) ,  ZO 

NAMELIST/DATA/ IPROJ .POLON. IBEAM, TITLE .LIMIT, LATMX , 

1  LATMN , LONMN , LONMX , POL AT , SCLLAT , 3CLL0N 

DECLARE  TAPE  FOR  CONTINENTAL  BACKGROUND  DATA  FILE 

OPEN (UN  I T =8 , FGRM=  UNFORMATTED') 

CALL  PRFLOT ( 3H300 , 4HPL0 T  > 

ID  <  1  >  *  1 
ID (2)  =  2 
ID  <  3)  =  -1 
ID  <  4)  =  0 
ID (5)  -  1 
ID (6)  -  2 
ID (7)  -  -1 

ZO  =  <0.1009,0.1076) 

WRITE  <6, 1 )  ZO 

1  FORMAT <2X,  LAT -LONG <RAD>  -  ,2E20.8> 

PI  *  3.1415926536 
S  =  1. 131463449 
A I  *  98. 1 
LTYPEl  *  2 
I BM  INI  ■=  21 
ISYMB1  *  7 
I FREQ 1  «  1 

Zl  *  < 180. /PI) *AIMA6(Z0> 

Z2  *  < 180. /PI ) *REAL  <  Z0> 

WRITE (12)  LTYPEl ,IBMIN1 , ISYMBl , I FREQ 1 , 1  ,Z1  ,Z2 
A ( 1 >  =  0.0047124 
A ( 2)  *  0.006109 
A  <3)  *  0.006807 
A  < 4 )  »  0.007697 
A <5)  *  0.008378 
A (6)  *  0.008901 
A (7)  ■  0.009599 
A (8)  3  0.010123 
A <9)  *  0.010472 
A ( 10)  -  0.011519 
A  < 1 1 )  *  0.012217 
A  ( 12)  x  0.012915 
A < 13)  -  0.013614 
A < 14)  -  0.014137 
A  (13)  -  0.014833 
A <  16)  *  0.015339 
A  <  17)  =  0.015708 
A  ( 18)  -  0.016381 
A  < 19)  -  0.017453 
A (20)  «  0.018850 
A <21 )  «  0.0202438 
A (22)  *  0.02094395 
A (23)  -  0.021293017 
A  (24)  x  0.021642 
A ( 25)  *  0.0219911 
A (26)  ■  0.0223402 
A (27)  -  0.022689 
A  <  28)  x  0.02740 
A (29)  x  0.0366319 


•• 

M  V  ^  J  /  Jd 

710  C 

SIGMA  =  (PI/4.)  +  0.9*0.01745329 

720 

SIGMA  *  (PI  /4.  ) 

730 

WRI TE  <6.2)  AI ,S,A(2) .SIGMA 

740 

2  FORMAT ( 2 X,  INCLINATION  OF  ORBIT iDEG)  = 

F 1 0 . 3 , / ,2X 

730 

1  SATELLITE  ALTITUDE  IN  EARTH  RADII  UNITS 

-  ,  F  1 0 .  3 

760 

2  '3-DB  BEAMWIDTH (RAD)  *  ' , F10. 6 , / , 2X , 

770 

3  ANTENNA  ROINTING  ANGLE (RAD)  =  .F10.3) 

700 

RO  -  (PI / 180. ) *AI 

790 

R1  =  COS(RO) 

000 

R2  »  S IN (R0> 

810 

R3C0S (SIGMA) 

820 

R4  -  SIN(SIGMA) 

830 

DBETA  *  S*R3  -  SORT ( 1 .  -  (S**2> * (R4**2) ) 

840 

WR I TE (6,3)  DBETA 

850 

3  FORMAT (2X,  DISTANCE  OF  SATELLITE  TO  SCAN 

POINT  =  ,1 

860 

LTYPE  *  1 

870 

IBM IN  =  17 

880 

ISYMB  =  0 

890 

IFREQ  *  1 

900 

SCLLAT  =  1. 

910 

SCLLON  *  1. 

920  C 

930  C 

TIME  LOOP 

940  C 

950 

DO  10  1=1, 200 

960 

II  =  1-1 

970  C 

980  C 

DECLARE  TAPE  FOR  LATITUDE  AND  LONGITUED  POINTS 

990  C 

1000 

OPEN ( UN I T= 1 2 , FORM-  UNFORMATTED  ) 

1010 

IF ( I . GT . 23)  GO  TO  70 

1020 

FI  =  II  ♦  265. 

1030 

GO  TO  80 

1040 

70  IF ( I . GT. 30)  GO  TO  71 

1050 

FI  -  MOD (II, 23 )  +  2165. 

1060 

GO  TO  80 

1070 

71  IF ( I . GT . 75)  GO  TO  72 

1080 

FI  -  MOD (II, 50 >  +  4065. 

1090 

GO  10  80 

1100 

72  IF ( I . GT. 100)  GO  TO  73 

1110 

FI  -  MOD (I I, 73)  ♦  3963. 

1120 

GO  TO  80 

1130 

73  IF ( I, GT. 123)  GO  TO  74 

1 1 40 

FI  =  MOD (I I, 100)  +  7865. 

1150 

GO  TO  80 

1160 

74  IF ( I . GT . 150)  GO  TO  75 

1170 

FI  -  MOD (11,1 23)  +  9765. 

1180 

GO  TO  80 

1190 

73  IF ( I . GT . 1 75)  GO  TO  76 

1200 

FI  -  MOD (11,1 30 )  *  11663. 

1210 

GO  TO  80 

1220 

76  FI  -  MOD (II, 173)  +  13565. 

1230 

80  CONTINUE 

1240 

R9  -  0.0009131 73*0. 001 *FI 

1 230 

RIO  -  COS (R9) 

1260 

Rll  -  S I N ( R9 i 

1270 

R12  -  Rll *R2 

1280 

R13  =  ACOS (R12) 

1290 

R14  -  SIN (R13) 

1300 

R13  -  R11*R1/R14 

1310 

R16  -  AS  IN (R15) 

1320 

R17  -  (PI/2. )  -  R13 

1330 

R18  -  ( 180. /PI ) *R17 

1340 

R19  -  ( 180. /PI >*R16 

1330 

WRITE (6,4)  R18.R19 
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1 3o0  4  r-ukfW  t  C»,  Lm 111  UbE  OF  ShI  ELl  11 1 ( DEG )  =  ,  EL'U .  fa ,  /  ,  _’A  , 

1370  1  LONGITUDE  OF  SATELLITE  DEGREES  WEST  =  .E20.8) 

1380  BETA  =  0.003306940*FI 

1390  R20  =  COS (BETA) 

1400  R21  =  SIN (BETA) 

1410  C 

1420  C  BEAMWIDTH  LQOP 

1430  C” 

1440  DO  40  L*1 , 30 

1450  R3  *  A (L) 

1460  R6  =  1.  ♦  R5**2 

1470  R7  ■  S**2  -  1. 

1480  R8  *  R6*R7 

1490  LASTJ  »  0 

1500  JJ  =  0 

1510  C 

1520  C  FOOTPRINT  GENERATOR  LOOP 

1530  C 

1540  DO  20  J»1 , 181 

1550  FJ  ««  2.  *J 

1560  OMEGA  =  (PI/180. )*FJ 

1570  R22  *  COS (OMEGA) 

1580  R23  *  SIN (OMEGA) 

1590  R24  *  -R3  ♦  R3*R22*R4 

1600  R25  *  -S*R24 

1610  R26  »  R25*R25  -  R8 

1620  IF (R26.GT.0. >  GO  TO  30 

1630  LASTJ  -  J 

1640  GO  TO  20 

1650  30  JJ  “  J  -  LASTJ 

1660  T  -  (R25  -  SORT (R26> > /R6 

1670  R27  *  S*Rl 1*R2 

1680  R28  *  -R3*R11*R2 

1690  R29  *  -R4*R21*R1 

1700  R30  *  R4*R20*R 1 0*R2 

1710  R31  *  R 1 1 *R2/R3 

1720  R32  -  R20*R1 

1730  R33  *  R21*R2*R10 

1740  R34  *  R3*R10 

1750  R35  *  R4*Rli*R20 

1760  R36  »  -R10/R3 

1770  R37  *  S*R21*Rl 1 

1780  R38  *  R1 1*R1 

1790  R39  »  -R3*R1 1*R1 

1800  R40  «  R4*R21*R2 

1810  R41  *  R4*R20*R10*R1 

1820  R42  -  R11*R1/R3 

1830  R43  «  R20*R2 

1840  R44  -  -R21*R1*R10 

1850  COSTHET  ■  R27  ♦  T*(R28  ♦  R29  ♦  R30  ♦  (R5*R22*R3/R4) * (R28  ♦ 

1860  1  R29  ♦  R30  ♦  R31)  -  <R5*R23) * (R32  +  R33) ) 

1870  STHCPH  -  S*R 10  -  T*<R34  ♦  R35  +  (R5*R22*R3/R4 > * (R34  ♦  R35  + 

1880  1  R36)  -  <  R5*R23 ) *R37 ) 

1890  THSF'H  =  R3B  ♦  T*(R39  +  R40  *  R41  ♦  (R5*R22*R3/R4>  *  (R39  ♦ 

1900  1  R40  ♦  R41  +  R42)  +  (R5*R23> * (R43  +  R44)> 

1910  R45  -  ACOS (COSTHET) 

1920  R46  ■  ATAN2(STHSPH, STHCPH) 

1930  R47  -  (PI/2. )  -  R45 

1940  W(J)  ■*  R46  +  (0. ,1.)*R47 

1950  IF ( I . GT. 84 )  GO  TO  20 

1960  UFQOT ( JJ , I )  =  R47/ (PI/180. ) 

1970  VFOOT ( J J , I )  ■  R46/ (PI/180. > 

1980  20  CONTINUE 

1990  IF( (L.EQ. 1. AND. I .LE.84) .OR. (L.EQ.9. AND. I.LE.B4) )  GO  TO  21 

2000  GO  TO  22 

2010  21  WRITE  (12)  LTYF'E  ,  IBM  IN,  ISYMB,  I  FREQ, 


MT^  %•  -P  wTimMV.W 


* 

JU20 

1  JJ  ,  (UFUOT  »►  ,  l  >  ,VHJUT  <K,  1  / 

20  ZO 

22  CONTINUE 

> 

2040  C 

2050  C 

CHECK  IF  Z0  IS  INSIDE 

- 

2060  C 

k" 

2070 

I SUM  -  0 

| 

2080 

DO  11  K-1,181 

! 

2090 

KK  -  K  +  1 

2100 

IF  <KK. EQ. 182)  KK-1 

!j 

2110 

KKK  *  INQUAD (W(KK>  -  Z0>  - 

2120 

JD  ■  ID <4  ♦  KKK) 

2130 

IF  ( JD. EQ. 2)  GO  TO  12 

2140 

11  1SUM  ■  I SUM  +  JD 

kl 

2150 

IF  < I SUM. EQ. 0)  GO  TO  12 

■ 

2160 

WR I TE  <6,13)  FI,L,A(L> 

2170 

13  FORMAT <2X,  TIME (MILLISEC) 

2180 

1  (L  *  ,15,')  *  ,F10. 7) 

1*. 

2190 

60  TO  10 

2200 

12  CONTINUE 

2210 

40  CONTINUE 

2220 

10  CONTINUE 

A 

2230 

CLOSE (12) 

1 

2240 

OPEN (12, FORM* ’ UNFORMATTED ' 

»  ' 

2250 

REWIND  12 

2260 

CALL  WORLDS (5,6) 

2270 

CALL  ENPLOT 

•» 

2280 

CLOSE (12) 

2290 

CLOSE (8) 

2300 

STOP 

| 

2310 

END 

1 

2320 

FUNCTION  INQUAD (ZZ> 

• 

2330 

COMPLEX  Z,ZZ 

► 

2340 

Z=ZZ 

V 

• 

2350 

IF (REAL ( Z ) )  1,2,3 

2360 

1  I F ( A I MAG ( Z ) )  11,11,12 

■ 

2370 

11  INQUAD  ■  3 

k 

2380 

RETURN 

V 

2390 

12  INQUAD  *  2 

.* 

2400 

RETURN 

j> 

2410 

2  I F ( A I MAG ( Z ) )  21,22,12 

2420 

21  INQUAD  -  4 

2430 

RETURN 

> 

2440 

22  INQUAD  =  0 

V 

l 

2450 

RETURN 

1 

2460 

3  IF (AIMAG ( Z ) )  21,31,31 

1 

2470 

31  INQUAD  *  1 

F 

2480 

RETURN 

K 

2490 

END 

**  •  »V*  »% 


list 


100  PROGRAM  F00T2  < INPUT , TAPES'  INPUT , OUTPUT , TAPE6=0UTF'UT ) 

110  DIMENSION  UFOOT (361 > , VFOQT (361 ) 

120  NAMELIST /DATA/ I PRO J ,  F'QLON, I  BEAM, T I TLE , L IMI T ,LATMX , LATMN, 

130  1  LONMN , LONMX , POL AT , SCLLAT , SSCLON 

140  OPEN (UNIT*8,F0RM= 'UNFORMATTED ' ) 

ISO  CALL  F'RPLOT  <  3H300 , 4HPL0T ) 

160  BEAMX  *  0.03 

170  BEAMY  =  0.04 

ISO  ROT  *=  0. 

190  PI  =  3. 1413926536 

200  S  -  6.619192 

210  SATLON  «  100. 

220  A I ML AT  -  40. 

230  AIMLON  -  100. 

240  LTYPE  «  1 

230  IBMIN  »  17 

260  ISYMB  »  O 

270  I  FREQ  -  1 

280  WRITE <6, 2>  BEAMX , BEAMY , ROT , SATLON, AIMLAT , AIMLON 

290  2  FORMAT <2X,  BEAMX <RAD>  -  ' ,F10. 3,2X ,  BEAMY (RAD)  «  ',F10.5,2X, 

300  1  'ROTATION (RAD)  *  ' ,F 10. 3 , / , 2X , ' SATLON (DEG)  ■= 

310  2  F 10.  3 , 2X  ,  '  AIMLAT  (DEG)  *  '  ,F  10.  3 , 2X  ,  '  AIMLON  ( DEG)  •«  ', 

320  3  F 10. 3) 

330  C 

340  C  PLOT  FOOTPRINT 

330  C 

360  RATIO  -  PI/ 180. 

370  FI  «  COS ( A I MLON*RAT 10) *COS ( A I ML AT *RAT 1 0 >  -  S*COS (SATLON* 

380  1  RATIO) 

390  F2  -  F 1*SIN ( AIMLAT*RATIO> 

400  F3  »  SIN ( AIMLON*RAT 10) *COS (AIMLAT *RATIO>  -  S*SIN (SATLON* 

410  1  RATIO) 

420  F4  »  F3*SIN  <AIMLAT*RATIO) 

430  A2  ■  S**2  -  1. 

440  D1  =  SORT ( 1 .  +  S**2  -  2. *S*COS (AIMLAT *RATIO) *COS (RATIO* ( 

430  1  AIMLON  -  SATLON) ) ) 

460  D2  ■*  SQRT (Dl**2  -  SIN(AIMLAT*RATIO) **2) 

470  LASTJ  -  0 

480  DO  20  J= 1,361 

490  FJ  -  J 

300  OMEGA  «  FJ*RATIO 

310  F3  »  BEAMX*COS (OMEGA) *COS (ROT)  -  BEAMY*SIN (OMEGA) *SIN (ROT) 

320  F6  «  BEAMX *COS (OMEGA) *S IN (ROT)  +  BEAMY *S IN (OMEGA) *COS (ROT) 

330  AO  -  1.  +  ( BEAM X*COS( OMEGA)  >**2  (BEAMY*SIN  (OMEGA)  )  **2 

340  A1  -  2.*S*(C0S(SATL0N*RATI0)*( (Fl/Dl)  -  (F2*F5/ ( D1 *D2) )  + 

350  1  (F3*F6/D2> >  ♦  SIN (SATLON*RAT 1 0) * ( (F3/D1)  -  (F4*F5/ (Dl*D2) ) 

360  2  -  (F1*F6/D2> ) ) 

570  F7  *  A 1**2  -  4. *A0*A2 

380  IF (F7. GT. 0. )  GO  TO  6 

390  LASTJ  =  J 

600  GO  TO  20 

610  6  JJ  -  J  -  LASTJ 

620  T  -  ABS (-A1  -  SQRT (F7) ) / (2. *A0> 

630  F 10  ■  T/Dl 

640  PHI  I  ■  ASIN(  (F 10)  *  (SI  N  ( AIMLAT*RAT  10  >  «•  F3*D2>  ) 

630  FB  -  S*SIN (SATLON*RAT ID)  +  T*((F3/D1)  -  (F4*F5/ ( D1 *D2 ) > 

660  1  -  (F1*F6/D2> ) 

670  F9  -  S*COS ( SATLON*RAT I 0 )  +  T* ( (Fl/Dl)  -  (F2*F5/ <D1 *D2) ) 

680  1  <F3*F6/D2)> 

690  ALAMDI  ■  ATAN2 (FB,F9) 

700  UFOOT(JJ)  *  F'HI  I  /RATIO 
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AlAMDI/RATIO 


/ iu  vFOOr  < jj>  = 

720  20  CONTINUE 

,30  WRITE (12)  LTYPE, IBMIN, IS/MB, IFREP 

740  1  JJ, (UFOOT(k) ,VFOOT(K> ,K=1,JJ> 

750  CLOSE (12) 

760  OPEN  <12, FORM= ' UNFORMATTED ' ) 

770  REWIJ2 

780  CALL  WORLDS (5,6) 

790  CALL  ENPLOT 

800  CLOSE (12) 

810  CLOSE (8) 

820  STOP 

830  END 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  "architect-engineer"  for 
national  security  projects,  specializing  in  advanced  military  space  systems. 
Providing  research  support,  the  corporation's  Laboratory  Operations  conducts 
experimental  and  theoretical  investigations  that  focus  on  the  application  of 
scientific  and  technical  advances  to  such  systems.  Vital  to  the  success  of 
these  investigations  is  the  technical  staff's  wide-ranging  expertise  and  its 
ability  to  stay  current  with  new  developments.  This  expertise  is  enhanced  by 
a  research  program  aimed  at  dealing  with  the  many  problems  associated  with 
rapidly  evolving  space  systems.  Contributing  their  capabilities  to  the 
research  effort  are  these  individual  laboratories: 

Aerophysics  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat 
transfer  and  flight  dynamics;  chemical  and  electric  propulsion,  propellant 
chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
spacecraft  structural  mechanics,  contamination,  thermal  and  structural 
control;  high  temperature  thermomechanics ,  gas  kinetics  and  radiation;  cw  and 
pulsed  chemical  and  excimer  laser  development  including  chemical  kinetics, 
spectroscopy ,  optical  resonators,  beam  control,  atmospheric  propagation,  laser 
effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions, 
atmospheric  optics,  Light  scattering,  state-specific  chemical  reactions  and 
radiative  signatures  of  missile  plumes,  sensor  out -of -f ield-of -view  rejection, 
applied  laser  spectroscopy,  laser  chemistry,  laser  optoelectronics,  solar  cell 
physics,  battery  electrochemistry,  space  vacuum  and  radiation  effects  on 
materials,  lubrication  and  surface  phenomena,  thermionic  emission,  photo¬ 
sensitive  materials  and  detectors,  atomic  frequency  standards,  and 
environmental  chemistry. 

Computer  Science  Laboratory:  Program  verification,  program  translation, 
performance-sensitive  system  design,  distributed  architectures  for  spaceborne 
computers,  fault -tolerant  computer  systems,  artificial  intelligence,  micro¬ 
electronics  applications,  communication  protocols,  and  computer  security. 

Electronics  Research  Laboratory:  Microelectronics,  solid-state  device 
physics,  compound  semiconductors,  radiation  hardening;  electro-optics,  quantum 
electronics,  solid-state  lasers,  optical  propagation  and  communications; 
microwave  semiconductor  devices,  raicrowave/raill imeter  wave  measurements, 
diagnostics  and  radioraetry,  microwave/millimeter  wave  thermionic  devices; 
atomic  time  and  frequency  standards;  antennas,  rf  systems,  electromagnetic 
propagation  phenomena,  space  communication  systems. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metals, 
alloys,  ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon;  non¬ 
destructive  evaluation,  component  failure  analysis  and  reliability;  fracture 
mechanics  and  stress  corrosion;  analysis  and  evaluation  of  materials  at 
cryogenic  and  elevated  temperatures  as  well  as  in  space  and  enemy-induced 
environments. 

Space  Sciences  Laboratory:  Magnetospheric,  auroral  and  cosmic  ray 
physics,  wave-particle  interactions,  magnetospheric  plasma  waves;  atmospheric 
and  ionospheric  physics,  density  and  composition  of  the  upper  atmosphere, 
remote  sensing  using  atmospheric  radiation;  solar  physics,  infrared  astronomy. 
Infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  ionosphere  and  magnetosphere; 
effects  of  electromagnetic  and  particulate  radiations  on  space  systems;  space 
instrumentat ion. 
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